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ABSTRACT 

In this letter we consider the evolution of a planetary system around a star inside a 
wide binary. We simulate numerically the evolution of the planetary orbits for both 
co-planar and highly- inclined systems. We find that the Kozai mechanism operates in 
the latter case. This produces a highly eccentric outer planet whose orbit crosses those 
of some of the inner planets. Strong planet-planet interactions then follow resulting in 
the ejection of one or more planets. We note that planetary systems resembling our 
solar system, formed around single stars in stellar clusters may exchange into binaries 
and thus will be vulnerable to planet stripping. This process will reduce the number 
of solar-system like planetary systems, and may produce at least some of the observed 
extra-solar planets. 

Key words: Celestial mechanics, Stellar dynamics; Binaries: general, Planetary sys- 
tems 



1 INTRODUCTION 

Of the 200 know n extrasolar planets (see for example 
iButler et af] 120061 ) around 40 l|Desidera fc Barbieril |2007h 
are in known binary systems, while the fraction with an 
unseen companion (e g. brown dwarf ) may also be signifi - 
cant (see for example Burgasser 20Qfj; iLuhman et alJ l2007). 
A key question is how the presence of the companion af- 
fects the orbits of the planets and the stability of the sys- 
tem. Various effects must be considered; if the systems 
consists of one planet inside a tight binary, a close en- 
counter between the planet and the companion can lead 
to the ejection of the planet (IH olman fc W icacrt 199!l). In 
wider binaries, secular effects, which builds up ove r thou- 
sands of orbits are important. IHolman fc WiegertJ (|l999l ) 
derived limits on the values which the semi-major axis and 
eccentricity of a companion star can have, without caus- 
ing large effec t s on t he orbits of the planets. Furthermore, 
M arzari et al l l|2005l ) analyzed the importance of planet- 
planet interactions in multi-planet systems with a co-planar 
companion star. If the companion star is highly inclined to 
the orbital plane of the planets, an effect called the Kozai 
Mechanism l|Kozailll962l ) becomes important. Its importance 
for single-planet sys t ems h as been analyzed in for exam- 
ple iTakeda fc Rasiol i|2005l ). Given that the solar system 
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is stable, it is possible to derive limits on the properties 
(mass and distance) that a hypothetical companion star 
to the sun could have, without causing severe damage to 
the s olar system (|lnnanen et alj 1 19971 ; IHolman fc Wiegertl 
1999). In this letter we extend these studies, consider- 
ing in more detail the evolution of planetary systems re- 
sembling our own solar system, when placed inside a bi- 
nary. Formation of planetary systems in binaries is possibl e 
and has been studied in for example Barbieri et al. (l2002t) ; 
iMarzari fc Scholll ll200o1);lThebault et all ( |2004 . l2006h : lBos3 



1 2006)7 rQuintana fc Lissauerl ( 20061 ). In this letter however, 
we only consider planet formation around single stars. The 
key idea here is to have such single stars exchange into bina- 
ries in young stellar clusters. In such three-body encounters, 
it is possible that the planetary system will be destroyed 
immediately, but a significant fraction will be unaffected, 
assuming that th e final binary is wide (see for example 
I Davies et al .11 19941 ). If this star had a planetary system this 
would now be randomly oriented with respect to the com- 
panion star. It is the evolution of these system that we study 
in this letter. In section [2] we discuss the integrator used to 
perform the numerical simulations presented in this letter. 
Then we give an overview of the Kozai mechanism in sec- 
tion [3] We present the results of our numerical simulations 
in section r?] which show that planet-planet interactions can 
cause the disruption of planetary systems if inside an in- 
clined binary. In section [5] we discuss the implications for 
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extra-solar planets and the vulnerability of planetary sys- 
tems resembling our solar system. 



2 SIMULATIONS - THE MERCURY CODE 

The numerical simulations in this paper where done using 
a modified version of the publicly available mercury code, 
which we have specifically adapted for integrating planetary 
systems in binary sta r s. Th e original mercury code is de- 
scribed in IChambersI (1 19991), and the modifications in our 
version are described in Chambers et al. 1 (I2002D 

The original mercury package (|Chamberslll999l ) can 
accurately integrate planetary systems, including the effects 
of close encounters between planets, very fast. The high per- 
formance is achieved because it intrinsically conserve the 
angular and linear momentum of the system. Integrators 
which have this feature, are said to be symplectic. Another 
property of them is that they do not show any long-term 
build up of energy error, as conventional integrators do, al- 
though they do exhibit high frequency energy oscillations. 
The original mercury package could not accurately simulate 
a planetary system which was inside a binary. A method for 
adapting the me rcury package, to allow for this, was how- 
ever presented in lChambers et all |2002l ) and we have here 
implemented those changes. 



3 THE KOZAI MECHANISM 



The Kozai mechanism was first discussed in iKozail (|l962h . 
He analyzed how the orbits of inclined asteroids in the so- 
lar system were influenced by Jupiter. In this case, the 
mass of the asteroid is much smaller than the mass of 
Jupiter. This is also true if we instead take the example of 
a Jupiter size extra-solar planet and a companion star and 
thus the same analysis hol ds for extra-solar pl anets within 
wide and inclined binaries ( Innanen et al. 1 119971 ) (though see 
jCarruba et all |2002| . 120031 ) for a corrected version of the 
Kozai equations). The Kozai mechanism causes the eccen- 
tricity of the planet to vary periodically, if the orbits of the 
planet and the companion star are sufficiently inclined. It 
has been used to ex plain the high eccentricity of the planet 
around 16 Cygni B (jHolman et al.lll997l ; iMazeh et allll997l ) 
and also for more gene ral studies of the eccentricities of plan- 
ets in binary systems l|Haghighipourl [20041 ; iTakeda fc Rasiol 
120051 ). Its importance for more general thre e-body systems 
has also been discuss ed in lFord et all (2000). Furthermore. 
ICarruba et al.l (|2004l ) gives a detailed discussion and analy- 
sis of the Kozai mechanism and use it to explain the orbits 
of irregular satellites in the solar system. 

Returning to the application of the Kozai mechanism in 
this letter, one observes that it is only significant if the ini- 
tial inclination, io, between the orbits of the planet and the 
companion star is greater than 39.23°(sin~ 1 (^/2/5)), which 
hereafter will be referred to as the critical angle, i cr it. If this 
is the case, the inclination of the planet will oscillate between 
icrit a nd ip and on e also observes that y/1 — e 2 cos (i) is con- 
stant l|Kozailll962l ). It can thus be shown, that the maximum 
eccen tricity, e max , of the planets orbit is (jlnnanen et alj 
Il997h : 




Figure 1. This plot shows the eccentricity evolution of a Neptune 
mass planet, with the same semi-major axis and eccentricity as 
Neptune (a = 30.08 AU and e = 0.0088), orbiting a one solar 
mass star within a binary system, where the properties of the 
companion star are: M c = 0.5Mq, a c = 300 AU, e c = 0.3 and 
i c = 60° . On the horizontal axis we plot the time in years and on 
the vertical axis the eccentricity of the planet. The triple system 
is formed at t = and the duration of the run spans 10 000 
orbital periods of Neptune. The oscillations observed are in good 
agreement with what we expect from the Kozai mechanism. 



where it is assumed that the the planets initial eccentricity, 
eo is close to zero. This is a lower limit on e max , since if eo 
is greater than 0, the maximum eccentricity will increase. 
Using equation [1] one can see that if for example io = 60°, 
the eccentricity will reach a maximum value of e max = 0.76. 

It is important to note that the effects of the 
Kozai mechanis m can be wa s hed o u t. The effects of non - 
spherical stars l|Soderhielm 1980l: Innanen et al. 1997 ). 
relativistic precession ( Holman et al.l 19971 ; Ford et all 
2000) and planet-pl anet interactions over long timescales 
(jlnnanen et al.lll997l ) can all decrease the eccentricity build- 
up. If the Kozai mechanism is to dominate over the planet- 
planet interactions, its timescale must be sufficiently short. 
This depends strongly on the semi-major axis of the com- 
panion star and the planet and thus, in order to see a large 
eccentricity build-up, the semi-major axis of the companion 
star must not be too large. 

To illustrate the Kozai mechanism, we ran simulations 
with the mercury code with a Neptune mass planet on 
Neptune's orbit (e p = 0.0088 and a p = 30.07 AU), orbiting 



e m ax = — 5/3 cos 2 (i ). 



(1) 



a one solar mass star, with a companion star with semi- 
major axis a c — 300 AU, eccentricity e c = 0.3, inclination 
io = 60° and mass M c = 0.5Mq. The parameters of the com- 
panion star are chosen to be representative for the known 
comp anions to planet-hosting s tars in the solar neighbor- 
hood (|Desidera fc Barbierill2007l ). 

We plot the evolution of the eccentricity of the planet in 
one of these runs in Fig.[T] As can be seen, its eccentricity e p 
oscillates, as predicted by the Kozai mechanism, with max- 
imum amplitude e max = 0.77. We can compare this to the 
theoretical value, given by equation[T] which gives that e m ax 
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should be equal to 0.76. Furthermore, one can see that the 
initial build-up of eccentricity is slow and then accelerates. 
It can be shown ()lnnanen et al]ll997l ) that the time to reach 
the first eccentricity maximum scales as T oc a,p 3 ^ 2 , where 
dp is the semi-major axis of the planet, under the assump- 
tion that the binary parameters and the planets eccentricity 
are kept constant. We can now compare the Kozai time- 
scales of for example Uranus and Neptune and find that it 
is about twice as long for Uranus as it is for Neptune. Thus, 
in a planetary system consisting of two or more planets, 
the eccentricity of the outer planet will reach its maximum 
when the eccentricity of the inner planets is still almost at 
its initial value. This is because the eccentricity stays at its 
original value for quite a long time until the growth starts, 
but once it does, it is very fast, as can be seen in Fig. [T] 

Thus, if an extra-solar planetary-system is within a 
highly inclined and not too wide binary, the Kozai mech- 
anism can lead to a very large increase in the eccentricity 
of the outer planet. When this occurs, the inner planets will 
still be almost unaffected by the companion star. If the in- 
clination is high enough, the eccentricity of the outer planet 
will be large enough for the orbits of the two outer planets to 
cross and hence strong planet-planet encounters are likely to 
occur. It is however very important to note, that when the 
eccentricity of the outer planet is at its maximum, the incli- 
nation between it and the companion star is at its minimum. 
This also means that the outer planet is inclined to the rest 
of the planets. This decreases the likelihood of planet-planet 
interactions. The outcome is that it can take several Kozai 
cycles before the planets have strong encounters (which de- 
pends on the phasing of the planets). Once they occur, such 
strong encounters can lead to large changes in both the semi- 
major axis and the eccentricity of the planets, or even the 
ejection of one or two of them, as we show in the following 
section. 



4 NUMERICAL SIMULATIONS - RESULTS 

As shown in the previous section, the Kozai mechanism can 
lead to strong planet-planet interactions. We investigate this 
here numerically, by taking the four gas giants of our own 
solar system around a one solar mass star and evolve the 
system within a binary. In Fig. [2] we plot the results of such 
a run, where the system has a companion star with m c = 
O.5M , a c = 300 AU, e c = 0.3 and i c = 60°. 

As can be seen, there is an almost immediate build-up 
in the eccentricity of the outer planet, and after only about 
3 million years the two outermost planets have been ejected, 
due to strong planet-planet interactions. By looking at the 
history of the close encounters that the two outer planets 
in the run in Fig. [2] had, it is possible to better understand 
what happens. 

Both of them have a series of close encounters with the 
other planets and themselves, which coincide with their pe- 
riods of high eccentricity over the first million years of the 
run. After this, there is a large change in the semi-major axis 
of the outermost planet, which leads to even more strong en- 
counters between all the planets. Both the outer planets are 
then ejected shortly after each other, at around 2.5 x 10 6 
years, due to strong encounters with the innermost planet. 
One can from the close encounter history clearly see that 
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Figure 2. This plot shows the evolution of the the four gas giants 
in our solar system, placed within a binary system. The properties 
of the companion star are: a c = 300 AU, e c = 0.3, m c = 0.5Mq 
and i c = 60°. The orbital elements of the companion star does not 
change during the duration of the simulation. On the horizontal 
axis we have plotted the time in years and on the vertical axis we 
have in the lower figure plotted the semi-major axis in AU and in 
the upper figure the eccentricity of the planets. 

all the sudden changes in the orbital elements of the planets 
seen in Fig. [2] are solely due to close encounters with other 
planets. 

It could be possible that the results seen in Fig. [5] are 
due to the initial phasing of the planets and that planetary 
systems normally are less sensitive. To test this, we per- 
formed 10 different realizations of the run in Fig. where 
we varied the initial phasing of the planets. The results of 
these runs were similar to those of the run shown in Fig. [2] 
In all of the 10 runs performed, one or two outer planets in 
the system were ejected within 6 x 10 6 years. It is important 
to note, that in all the 10 runs, the companion star suffers 
no close encounter with any of the planets, hence the ejec- 
tion of planets are not due to strong interactions between 
the planets and the star. 

To further test our hypothesis that the instabilities seen 
in Fig. [2] are caused by the Kozai mechanism, through the 
increased eccentricity of the outer planet and the following 
planet-planet interactions, we then performed another set 
of 10 runs. These had the same initial conditions for the 
companion star as in Fig. [5] but with the initial inclination 
between the planets and the companion star equal to 10°. 
Hence, the initial inclination between the orbital plane of the 
binary and that of the planets is much less than the critical 
angle (39.23°), and so the Kozai mechanism should not be 
active. In all the 10 runs (which lasted for 10 8 years), there 
is no build-up of eccentricity for the outer planet and thus 
we see no strong planet-planet encounters, and no ejections, 
as is expected from the Kozai theory. 

If it is only planet-planet interactions that cause the 
ejection of one or more planets in the runs we performed, the 
outcome should be similar if we remove the companion star 
once the eccentricity of the outer planet is high. To test this, 
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we removed the companion star at a time when the eccentric- 
ity is at its maximum in the 10 highly-inclined runs discussed 
above. As expected, we still see planet ejections within the 
same time-scales, although the strong planet-planet encoun- 
ters that lead to the ejections are not the same, which is also 
expected. 

We have also performed simulations, where we place the 
four giants in a very wide binary to see whether the Kozai 
mechanism could be observed. This should not be the case, 
since the Kozai timescale becomes very long in this case, and 
its effect should be suppressed by the precessi on of the plan- 
ets d ue to planet-planet interactions (see also llnnanen et ah! 
1 19971 ). We again used the giant planets of the solar system 
as our test object, and this time placed a companion star 
at 1000 AU on a nearly circular orbit with an inclination of 
60°. The simulation showed that no secular evolution due 
to the effects of the companion star occurred during the 10 s 
years that we ran the simulation for. 



5 DISCUSSION 

As we have seen, a planetary system can be strongly affected 
by the presence of a companion star, even if the semi-major 
axis of the companions orbit is large. However, this is only 
true if the initial inclination between the orbital planes of 
the planets and the companion star is larger than the crit- 
ical angle, 39.23°. If we assume that all planetary systems 
that are in binaries also were formed in them, it would be 
unlikely to find systems with such a high mutual inclina- 
tion. This is because one would expect the rotational axis of 
the stars to be aligned and hence the disk out of which the 
planets form would be in the same plane as the companion 
star. However, if the planetary system is formed around a 
single star, which is later exchanged into a binary, the orien- 
tation of the planets orbital planes will be random. One can 
show, that for systems who have formed this way, 77% will 
have an initial inclination above 39.23° and hence be in the 
region where the Kozai mechanism can be important. To un- 
derstand how important the Kozai mechanism is for extra- 
solar planets, it is thus n ecessary to know how often stars are 
exchanged into binaries dPortegies Zwart fc McMillan! 120051 ; 
I Adams et all 120061 ; iPfahl fc Muterspaughl 120061 ). Stars are 
formed in groups. Thus, early on, stellar encounters will be 
frequent in such crowded environments. This means that it 
might be quite common for stars to be exchanged in and out 
of binaries. It is therefore interesting to study the fraction 
of single stars which have never been in a binary or even 
suffered a close encounter with another star. We term such 
stars si ngletons. The fraction of single stars has been cal- 
culated (.M almbe rg et al. 2007]) and the results show that a 
substantial fraction of the stars are not singletons. 

To understand how important such interactions be- 
tween stars are for planetary systems, we need a criterion for 
when they can cause significant damage to them. In the case 
that the star is exchanged into a binary for several 10 5 years, 
a companion star can cause several planets to be ejected, via 
planet-planet interactions induced by the Kozai mechanism. 
In the example of the solar system, Neptune's and Uranus 
orbits will cross if the inclination between the companion 
star and the plane of the planets is greater than 43°. This 
will be the case in 73% of the binary systems formed in 



stellar encounters. It is thus possible that planetary systems 
initially resembling our own solar system may more closely 
resemble the observed extra-solar planetary systems due to 
planet-planet scattering induced by the Kozai mechanism. 
This is true also for single stars with planetary systems, 
since these may have previously been in a binary and later 
exchanged out of it. 

If it is common for young planetary systems to be 
stripped of one or more planets while they are in young 
stellar clusters, we should be able to observe some of these 
planets. It is interesting to note that free-floating planetary 
mass object s has been observed in dense stellar clusters (see 
for example iLucas fc Roche] l2000l ). These objects might be 
young planets which have been ejected from its host sys- 
tem, but may also have bee n formed in other ways (eg. 
IWhitworth fc Zinneckerll2004 ) 

From the simulations which we have performed, we do 
not see the formation of so-called hot Jupiters. These are 
Jupiter mass planets found at very small distances (a few 
stellar radii) from the host star, on almost circular orbits. 
However, hot Jupiters may be formed from planets that 
have strong tidal encounters with the host star, which how- 
ever requires that the planets have very eccentric orbits 
(|Faber et al.l 1200*51 '). If the inclination between the planets 
and the companion star is near 90° the Kozai mechanism 
will, in the absence of relativistic effects, cause such high 
eccentricities, if the binary is not too wide. Thus, it might 
be possible to explain hot Jupiters in this way. However, 
observations imply that at least some highly eccentric plan- 
etary systems have avoided tidal circulation (|Wu fc Murray! 
2003). We will return to this problem in a subsequent paper. 

As we have seen above, no strong secular effects occur in 
our runs with planetary systems within a co-planar binary. 
This is true as long as the separation between the companion 
star and the central star is large. However, if the separation 
is smaller than 10 AU, one or more pl anets might be ejected 
from the system dMarzari et al.l [20051 ). Here, we do not see 
any significant build-up of the planets eccentricities in the 
low inclination case. However, we do see a small increase 
in the mutual inclinations between the planets and their 
respective eccentricities. 



6 SUMMARY 

We have simulated the evolution of planetary systems, re- 
sembling our own solar system, inside wide binaries. We find, 
that if the inclination between the orbital plane of the com- 
panion star and the orbital planes of the planets is larger 
than 39.23°, the Kozai mechanism will lead to an increase 
in the eccentricity of the outer planet, if the binary is not to 
wide. The increased eccentricity of the outer planet leads to 
strong planet-planet interactions in the system, which can 
lead to the ejection of one or more planets and also that the 
remaining planets are left on more eccentric orbits than be- 
fore. As stars are formed in groups, they will, during the first 
100 million years after their birth, be in a crowded environ- 
ment. This means that they can exchange into binaries in 
stellar encounters, with a random orientation of the orbital 
plane of the planets, with respect to the orbital plane of the 
companion star. Thus, a planetary system resembling our 
own solar system formed around a single star, is at risk of 
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exchanging into a binary and being stripped of one or more 
planets. 
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